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Abstract 

The lateral and longitudinal profiles of the hadronic showers de- 
tected by iron-scintillator tile hadron calorimeter with longitudinal 
tile configuration have been investigated. The results are based on 
100 GeV pion beam data. Due to the beam scan provided many 
different beam impact locations with cells it is succeeded to ob- 
tain detailed picture of transverse shower behavior. The underlying 
radial energy densities for four depths and for overall calorimeter 
have been reconstructed. The three-dimensional hadronic shower 
parametrisation have been suggested. 



1 Introduction 



Hadronic shower is a basic notion of hadron calorimetry. But despite that 
hadronic shower characteristics are studied for many years the exhaustive 
quantitative understanding of hadronic shower properties is not exist. The 
pubhshed data rule the energy deposition in calorimeter cells and 

therefore are related with specific cell dimensions and the acceptance of 
cells relative to shower axis. Furthermore, as to the transverse profiles they 
are as usual the energy depositions as a function of transverse coordinates, 
not a radius, and integrated over the other coordinate 0. Meanwhile for 
many purposes of experiments a very detailed simulation is not needed and 
a three dimensional parameterisation of hadron shower development is to 
become very important for fast simulation which significantly (up to 10^ 
times) to speed a detailed GEANT based simulation 0, 0, p. 

In this paper we report on the results of the experimental study of 
hadronic shower profiles detected by prototype of ATLAS tile hadron 
calorimeter This calorimeter has innovative concept of longitudinal 
segmentation of active and passive layers (see Fig. 0) and the measurement 
of hadron shower profiles therefore a special interest [^. This investiga- 
tion was performed on the basis of data from 100 GeV pion exposure of 
the prototype calorimeter at the CERN SPS at different z impact points 
in the range from —36 to 20 cm (2; scan) at incident angle O = 10° which 
were obtained in May 1995. 

Earlier some results related with lateral shower profiles for this calori- 
meter were obtained in 

2 The Calorimeter 

The prototype the ATLAS hadron tile calorimeter (Fig. 0) is composed of 
five modules. Each module spans 27r/64 in azimuthal angle, 100 cm in the 
z direction, 180 cm in the radial direction (about 9 interaction lengths A/ at 
?7 = or to about 80 effective radiation length Xq), and has a front face of 
100 X 20 cm^ The iron structure of each module consists of 57 repeated 
"periods" . Each period is 18 mm thick and consists of four layers. The first 
and third layers are formed by large trapezoidal steel plates (master plates) , 
5 mm thick and spanning the full radial dimension of the module. In the 
second and fourth layers, smaller trapezoidal steel plates (spacer plates) 
and scintillator tiles alternate along the radial direction. These layers 
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Figure 1: Principle of the tile hadronic calorimeter. 



consist of 18 different trapezoids of steel or scintillator, each spanning 
100 mm along x depending on their radial position. The spacer plates 
and scintillator tiles are 4 mm and 3 mm thick respectively. The iron to 
scintillator ratio is 4.67 : 1 by volume. The calorimeter thickness along 
X direction at incidence angle = 10° corresponds to 1.49 m of iron 
equivalent . 

Radially oriented WLS fibers collect light from the tiles at both of their 
open edges and bring it to photo-multipliers (PMTs) at the periphery of 
the calorimeter. Each PMT views a specific group of tiles, through the 
corresponding bundle of fibers. The calorimeter is radially segmented into 
four depth segments by grouping fibers from different tiles. As a result of 
each module is divided on 5 (along z) x A (along x) separate cells. 

The readout cells have lateral dimensions of 200 mm (along z) x 
(200 -T- 300) mm (along ?/, depending from a depth number) and longi- 
tudinal dimensions of 300, 400, 500, 600 mm for depths 1-4, corresponding 
to 1.5, 2, 2.5 and 3 A/ at r/ = 0. On the output we have for each event 200 
values of energies from PMT properly calibrated with pedestal 

subtracted. Here indexes i, j, k, 1 mean: i = 1, . . . , 5 is the column of cells 
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Muon wall 5 Module Detector 



Figure 2: Setup of the tile hadronic calorimeter prototype. 



number (tower), j = 1, . . . , 5 is the row (module) number, k = 1, ... ,4 is 
the depth number and / = 1, 2 is the PMT number. 

The calorimeter has been positioned on a scanning table, able to allow 
high precision movements along any direction. Upstream of the calorime- 
ter, a trigger counter telescope was installed, defining a beam spot of 2 cm 
diameter. Two delay-line wire chambers, each with {z,y) readout, allowed 
to reconstruct the impact point of beam particles on the calorimeter face 
to better than ± 1 mm . For the measurements of the hadronic shower 



longitudinal and lateral leakages backward (80x80 cm^) and side (40x115 
cm^) "muon walls" were placed behind and side the calorimeter modules 
i- 

Construction and performance of ATLAS iron-scintillator barrel had- 
ron prototype calorimeter is described elsewhere 0,0, fTH, [[T^ . 



3 The methods for extracting the underly- 
ing radial energy density 

The incident particle interacts with the material of calorimeter and a 
shower is initiated. In the following the used coordinate system will be 
based on the hadron shower direction. The axis of the shower, defined 
to be a track of a incident particle, forms the x axis with x = at the 
calorimeter front face. 

We measure the energy depositions in calorimetry cells. In ijk-cell of 
the calorimeter with volume Vijk and cell center coordinates {x, y, z) the 
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energy deposition Eijk is 

Eijk{x,y,z)= / f{x,y,z)dxdydz, (1) 

where f(x,y,z) is the three-dimensional hadron shower energy density func- 
tion. 

In the following we will consider the various marginal (integrated) den- 
sities longitudinal density fi{x) obtained after integrating f{x,y,z) 



over y and z coordinates 



Cymax pZmax 

fi{x) = / / f{x,y,z)dydz, (2) 



transverse densities f2{y,z) for four depths and overall calorimeter ob- 
tained after integrating f{x,y, z) over the various x ranges 

f2{y,z) = / f{x,y,z)dx, (3) 

J X\ 

transverse densities f{z) 

f{z) = / f{x,y,z)dxdy. (4) 

Cumulative function is 

F{z)= / f2{z,y)dzdy (5) 



oo J — oo 



and related with marginal density f{z) by 

/(.) ^ (6) 

Due to the azimuthal symmetry densities f2{y, z) are the function only 
a radius r = \J\y'^ + z^) from the shower axis, i.e. 

f2{y,z) = ^{r). (7) 

In such case 

Eijk =11 Hr)rdrd(l), (8) 
where (p is the azimuthal angle. 
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There are some methods for extracting of radial density $(r) from the 
measured distributions of energy depositions Eij^. 

One method is the unfolding $(r) from (§). This method was used in 
the analysis of data from the lead-scintillating fiber Spaghetti Calorimeter 
151 1 ■ Several analytic forms of $(r) were tried, but the simplest that 



describes the energy deposition in cells was a combination of an exponential 
and a Gaussian: 

$(r) = -ie ^1 + -^e ^^2> . (9) 

In order to determine the free parameters Oj, Ai in the expression a 
minimisation fit have been done. 

Another method is using the marginal density function f{z) and its 
connection with radial density $(r) | ]TB[ . 

^ f°° ^(r)rdr , , 

/(^)=2/, (10) 
— 



This method was used []TB[ for extracting of electron shower transverse pro- 
file on the basis of the data from GAMS'— 2000 electromagnetic calorimeter 
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Integral equation ([T0| ) can be reduced to the Abelian equation by re- 
placing of variables [|l^. We solved equation (|10D (see Appendix 1) and 
obtained 

1 d f°° f{z)dz^ 



^^^^ ndr^Jr^ y^;22-r2)' ^^^"^ 
The marginal density f{z) also may be determined by various ap- 



proaches. One approach is the using of cumulative function f{z) [jl6[ and 
differentiation of it according to (^. Another approach is to assume some 
form of marginal density f{z), to derive a formula for energy deposition in 
tower and to determine parameters of f{z) by fitting [|19 . 



If the three exponential distribution for parametrisation of f{z) is used 

i=l 

then for the energy deposition in a tower E{z), the cumulative function 
F[z), for radial density $(r) we obtain: 

/•z+h/2 

E{z)= f{z'-z)dz', (13) 

Jz-h/2 
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E{z) 
Eiz) 



f EL - e-^cMg)), for \z\ < ^, (14) 



f ELi«^A,s/i(2^)-e 



for \z\ > 



h 



(15) 



F{zr)= §ELl«^A,e'^^, forzr<0, (16) 

F(^r)= f + §ELia,A,(l-e-t), /or^,>0, (17) 

*«-i|:^Ao(£). (18) 

where z is the transverse coordinate of the tower center, Zr is the right edge 
transverse coordinate of the tower in the center, h is the front face size of 
towers along Z axis {h = 200 mm in our case), B = Ef=i (^iK, Eo, ai, Aj 
are free parameters. Eg is an energy normahsation factor, Kq is modified 
Bessel function. Using the condition Ei=i ctj = 1 we reduced the number 
of parameters to six. 

The radial containment of shower as a function of r is 

/(r)= r r$(r)rrfr^i0 = §^aA(l-fi^i(f)), (19) 
Jo JO ^ i=i Aj A, 

where Ki is modified Bessel function. 



4 Results 



Using the program TILEMON 30 runs with different impact points of 
incident particles in the range from —36 to 20 cm contained 320 K events 
have been analysed and the various energy deposition spectra have been 
obtained. 



4.1 Transverse behaviour of hadron showers 
4.1.1 Energy deposition in the towers 

Fig. 1^ shows the energy depositions into towers for 1 4 depths as a func- 
tion of z coordinate of center of tower. Fig. ^ shows the same for overall 
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Figure 3: The energy depositions in towers of 1 4 depths as a function 
of Z coordinate. Depth 1 — up left, depth 2 — up right, depth 3 — down 
left, depth 4 — down right. Curves are fits of equation (14) and (0) to 
the data. 
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Figure 4: The energy depositions in towers summed overall calorimeter 
depths as a function of Z coordinate. Curve is the result of fit by formula 
(14) and (p^Sl). Open circles are the Monte Carlo predictions (GEANT- 
FLUKA+GHEISHA) given in Table |. 
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calorimeter (the sum of histograms presented in Fig. |^). Due to the wide 
beam scan provided many different beam impact locations with cells and 
using information from all cells it is succeeded to obtain detailed picture of 
transverse shower behaviour in calorimeter. The obtained spread of trans- 
verse shower dimensions is more than 1000 mm. The energy depositions 
span a range of three orders of magnitude. It can be seen the fiat shoulder 
at 1 2 1 coordinate less than 100 mm. An immediate turnover occur as soon 
as 1 2; I reached the boundary of the cell. Such picture of transverse shower 



behaviour was observed in other calorimeters as well ||15[ , . 

We used these distributions in order to extract the underline marginal 
density f{z). By adding to statistical errors the electronic noise errors of 
27 MeV/cell the effective intercalibration error of 2% and uncer- 
tainties of 4%, arising from the nonzero entry angle of the incident beam 
into the calorimeter, we obtained a good description of these distributions. 

The solid curves in Fig. ^ and ^ are the results of fit with equation 



(|T5|). In compare with |T^, where the transverse profiles exists only 



for distances less than 250 mm, the our more extended profiles (up to 650 
mm) demand to introduce the third exponential. 

The parameters and Aj obtained by fitting are listed in Table |1]. As 
can be seen from Table |l| the slopes of exponentials, Aj, increase and the 
contribution of the first exponential, ai, decrease as shower develop. The 
obtained values of Ai and A2 for overall calorimeter agree well with ones 



for conventional iron-scintillators calorimeter amount to 18 ± 3 mm 
and 57 ± 4 mm, respectively. 



Table 1: The parameters and Aj obtained by fitting the transverse 
shower profiles for four depths and overall calorimeter. 





ai 


Ai, mm 


0-2 


A2, mm 


ag 


A3, mm 


1 


0.88 ±0.07 


17 ±2 


0.12 ±0.07 


48± 14 


0.004 ±0.002 


430 ± 240 


2 


0.79 ±0.06 


25 ±2 


0.20 ±0.06 


52 ±6 


0.014 ±0.006 


220 ± 40 


3 


0.69 ±0.03 


32 ±8 


0.28 ±0.03 


71 ±13 


0.029 ±0.005 


280 ± 30 


4 


0.41 ±0.05 


51 ± 10 


0.52 ±0.06 


73± 18 


0.07 ±0.03 


380 ± 140 


all 


0.78 ±0.08 


23 ± 1 


0.20 ±0.08 


58 ±4 


0.015 ±0.004 


290 ± 40 



The shower depth dependences of the parameters Oj and Aj are dis- 
played in Fig. |[ As can be seen they demonstrated a linear behaviour. 
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The curves are fits of linear equations ai = ai + I3i ■ x and Xi = 'ji + Si ■ x. 
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Figure 5: Tlie sliower deptli dependences of tlie parameters Oj and Aj. 

The values of parameters 7^ and 5i are presented in the Table |[ 

It is interesting that in for the low- density-fine-grained flash chamber 
calorimeter the linear behaviour of the slope exponential is also observed. 
In the same time for uranium-scintillator ZEUS calorimeter some non- 
linear behaviour of slope of halo component measured at 100 GeV has been 
demonstrated at interaction lengths more then 5 A in shower development 
21. 



Table 2: The values of parameters a^, 7^ and 5i. 





di 


l/K 




7i, mm 


5j, mm/XT, 


ai 


0.99 ±0.06 


-0.089 ±0.015 


Ai 


14 ±2 


6± 1 


a2 


0.04 ±0.06 


0.072 ±0.015 


A2 


42 ± 10 


6±4 


as 


-0.001 ±0.002 


0.008 ±0.002 


A3 


170 ± 80 


30 ±24 



Table ^presents the Monte-Carlo predictions for hadron shower energy 
depositions in the modules of TILECAL at centre coordinates of these 
modules obtained by using the different hadronic simulation packages in- 
terfaced with GEANT ||^ . The simulation energy cutoff for neutrons have 
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been set to 1 MeV. One can compare our data from Fig. ^ and Monte- 
Carlo predictions from Table ^ for the same z coordinates. Comparison 
shows worse agreement, with factor > 3, at distances of ±400 mm from 
shower kernel. 



Table 3: The average fractions of the energy deposited per module for 100 
GeV pions obtained by using the different hadronic simulation packages 
interfaced with GEANT (G). 



Z, mm 


-400 


-200 





200 


400 




E, % 


G-FLUKA+GHEISHA 


0.43 ± .02 


4.3 ± .1 


91 ± 1 


4.4 ± .1 


0.43 ± .02 


G-FLUKA+MICAP 


0.20 ± .01 


3.2 ± .1 


94 ± 1 


3.2 ± .1 


0.18 ± .01 


G-GHEISHA 


0.34 ± .01 


3.2 ± .1 


92 ±1 


4.5 ± .1 


0.38 ± .01 



4.1.2 Cumulative function 

Similar results were obtained from cumulative function distributions. Cu- 
mulative function F{z) was obtained as follows: 

F{z) = J:f\z), (20) 

k=l 

where F'^[z) is the cumulative function for fc-depth. For each event F^{z) 
is 

F\z)=Y.T.E:,k. (21) 

i=i j=i 

where imax = 1, . . . , 5 is the last number of tower in sum. 

Fig. H and ^ (left) present the cumulative functions F^{z) for four 
depths and for overall calorimeter. The curves are fits of equations (|l^) to 
the data. The results of the cumulative function fits are less reliable and 
in the following we will use the results from energy depositions in tower. 

The knowledge of the cumulative function allows us to determine di- 
rectly according to (^ the marginal density f{z) without additional as- 
sumption about its form. This is demonstrated in Fig. |^ (right) where the 
marginal density f{z) extracted by the numerical differentiation F[z) is 
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Figure 6: The cumulative functions -^(2;) for four depths. Depth 1 — up 
left, depth 2 — up right, depth 3 — down left, depth 4 — down right. 
Curves are fits of equations (16) and (0) to the data. 
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Figure 7: Left: The cumulative function f{z) for overall calorimeter. 
Curves are fits of equations (16) and ( [T^ ) to the data. Right: The marginal 
density f{z) for overall calorimeter. 

shown. Overlayed curve is a calculation of the expression (|12]) with pa- 
rameters listed in Table |I|. This curve well reproduces the such obtained 
marginal density. 

Thus, the marginal densities f{z) determined by three methods (by 
using the energy deposition spectrum, the cumulative function and the 
numerical differentiation of F{z) are in reasonable agreement. 

4.1.3 Radial hadron shower energy density 

Fig. 1^ and Fig. ^ (left) show the radial shower energy density functions, 
$(r), calculated by formula (p!8D with parameters from table |l]. The contri- 
butions of different terms are also shown. The data of SPACAL calorimeter 
calculated by formula (^ with Ai = 140 mm and A2 = 42.4 mm are shown 
for comparison on Fig. ^ (right). Since in this work the lateral profiles 
are presented in term of the measured picocoulombs the density $(r) was 
transformed into GeV energy scale by using the energy deposit constant 
oiApC/GeV !§. 

As can be seen at general reasonable agreement the curve of SPACAL 
lies systematically below the TILECAL one beyond 1.5 A^^-^. 
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Figure 8: The radial energy density, $(r), as a function of r for TILECAL 
for various depths: sohd hues — the energy densities $(r), dashed hues — 
the contribution from the first term, dash-dotted hnes — the contribution 
from the second term, dotted hnes — the contribution from the third term. 
Depth 1 — up left, depth 2 — up right, depth 3 — down left, depth 4 — 
down right. 
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Figure 9: Left: The radial energy density as a function of r (in units of 
A/) for TILECAL (solid line), the contribution to $(r) from the first term 
(dashed line), the contribution to $(r) from the second term (dash-dotted 
line), the contribution to $(r) from the third term (dotted fine). Right: 
The comparison of the radial energy densities as a function of r (in units 
of Xf) for TILECAL (sohd line) and SPACAL (dash-dotted hne). 
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4.1.4 Radial containment 
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Figure 10: Left: The containment of sliower /(r) (solid line) as a func- 
tion of radius for overall TILECAL calorimeter. Dash-dotted line — the 
contribution from first depth, dashed line — the contribution from second 
depth, thick dotted line — the contribution from third depth, thin dotted 
line — the contribution from fourth depth. Right: The radii of contained 
cylinders for given shower containment as a function of depths: black cir- 
cle - 90% of containment, black square - 95%, black triangle - 99%. The 
curves are drawn to guide the eye. 



One of the important questions concerns the shower transverse dimen- 
sions and its longitudinal development. 

The parameterisation of the radial density function, ^{r), have been 
integrated to yield the shower containment as a function of radius, I{r). 
Fig. |TD| (left) shows the transverse containment of the pion shower, /(r), 
as a function of r for four depths and overall calorimeter. 

In Table ^ and Fig. |l^ (right) the radii of cylinders for given shower 
containment (90%, 95%, 99%) extracted from Fig. |10| (left) as a function 
of depth are shown. Solid lines are the linear fit to the data. As can 
be seen these containment radii increase linearly with the depth. The lin- 
ear increase of 95depth is also observed for Fe-scintillator calorimeter at 50 
and 140 GeV For overall TILECAL calorimeter the 99% containment 
radius is equal to 1.7 ±0.1 A^'^-^. In the last row of this Table the corre- 
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Table 4: The radii of contained cylinders for given shower containment for 
TILECAL (for various depths and overall) and SPACAL (overall) calorime- 
ters. 



Calorimeter 


X in 


r, in X'J^ 


90 % 


95 % 


99 % 


TILECAL 


0.59 


0.44 


0.64 


1.43 


1.97 


0.60 


0.88 


1.55 


3.74 


0.92 


1.27 


1.75 


5.92 


1.24 


1.55 


1.87 


overall 


0.72 


1.04 


1.67 


SPACAL 


overall 


0.86 


1.19 


1.72 



sponding data for SPACAL calorimeter for 80 GeV vr-mesons are given (a 
pions grid scan at an angle of 2° with respect to fiber direction). Note that 
in the case of Pb-scintillator calorimeter SPACAL, having the same X^J^ 
(see Table |^ in Appendix 2), the shower containment radii are similar to 
obtained for iron-scintillator calorimeter TILECAL. 

It is interesting to note that it is mistaken to consider as the measure 
of the transverse shower containment the one obtained from the marginal 
density f{z) or the energy depositions in strips as have been made in 
In this work for 99% containment have been obtained the value of 1.2 A/ 
and the conclusion have been made that "their result is consistent with 
the "rule of thumb" that a shower is contained within a cylinder of radius 
equal to the interaction length of a calorimeter material [Q". But it is 
showed by our and SPACAL measurements that the value of radius of 
99% contained shower cylinder amounts to about two interaction length. 
Reduced value of this radius obtained from f{z) is due to that according to 
(P^) represents the integrated function $(r). In our case if we will use as a 
measure of shower containment the half- width of integrated f{z) estimated 
from F{z) in Fig. |^ (left) we obtain the value of 300 mm or 1.2 X'^J^ that 
agrees with [|T| . 

4.2 Longitudinal Profile 

Here we are concerned with the differential deposition of energy /S.E / Ax as 
a function of x, the distance along the shower axis. In Table |^ the average 
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energy shower depositions in various depths, Eg, the normahsed to one 
interaction length X^^ energy depositions, AE/Ax, the lengths of depths, 
L, the effective iron lengths of depths, Lf,ff, the lengths of depths in units 
of A^^, the centers of depth intervals in units of A^*^, x, are given. In these 
calculations the value of X^^ = 207 mm has been used (see Appendix 3). 



Table 5: Average energy shower depositions in various depths. 



a; in A^^ 


Ax in A^*^ 


L, mm 


Lpe, iiim 


Eo, GeV 


AE/Ax, GeV 


0.59 


1.18 


300 


245 


25.5±0.3 


21.6±0.3 


1.97 


1.58 


400 


327 


43.6±0.2 


27.6±0.1 


3.74 


1.97 


500 


408 


22.4±0.1 


11.2±0.1 


5.92 


2.37 


600 


490 


8.5±0.5 


3.6±0.2 



In Fig. |TT| the our quantities AE/Ax (open circle) together with the 
data of ||25| (open triangles) and Monte Carlo predictions {GEANT — 
FLUKA + MICAP, diamonds) are shown. The agreement is observed. 
So, as to longitudinal energy deposition our calorimeter with longitudinal 
orientation of the scintillating tiles agrees with the one for a conventional 
iron-scintillator calorimeters. 



4.3 "Electromagnetic" fraction of a hadronic shower 



The lateral shower profile information may be used for determination of the 
average fraction of energy going into 7r° production in a hadronic shower, 
/tto |T^. Following |]TB[ we assume that the "electromagnetic" part of 
hadronic shower is the prominent central core, i.e. in our case the first term. 
The integrated contribution from this term is fj^o = (55 ± 3)%. This value 
may be compared with the one of 56% as simulated by FLUKA for our 
calorimeter P|, the one of 59% as simulated by GALOR for iron calorimeter 



[ 26| and the one of 54% as obtained by using the lead-scintillating fiber 
SPAGAL calorimeter 
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The observed tt° fraction, /^o, is related to the intrinsic actual fraction, 
by the relation p^, 121 



flT° 



e/h ■ 



{e/h-l)-K. + l- 



(22) 
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Figure 11: Left: Longitudinal profile as a .function of longitudinal co- 
ordinate X in units A^^. Open circles are our data, open triangles are 
conventional calorimeter data, diamonds are the Monte Carlo predictions. 
Right: The fractions of the "electromagnetic" parts of showers in various 
depths as a function of x. 



There are two analytic forms for the intrinsic tt° fraction suggested by 
Groom ||2^ 

R (m-l) 

Jlo = l- (t^) (23) 



and Wigmans 



E 



Ko = k- ln{-) 



E' 



(24) 



where E'^ = 1 GeV, m = 0.85, k = 0.11. We calculated the values of 
using a value ofe/h= 1.23 for our calorimeter at = 10° |T2[ and obtained 
the values of /,r° are equal to 55% and 56% for Groom and Wigmans 
parameterisations of f^o, respectively. Thus, our measured value of /,r° 
agrees well with the experimental one 
26i and with our calculations. 



151 and Monte Carlo calculation 



The determined contributions of parts for various depths give the possi- 
bility to obtain the fractions of "electromagnetic" and "hadronic" parts of 
hadronic showers in various stages of longitudinal development of showers. 
On the Fig. |TT| (right) shows the corresponding results. As can be seen as 
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shower developed the "electromagnetic" fraction decreases. This is natural 
since the shower energy exhausts and as a result tt° production decreases. 
The best fit to these data is fem = (76 ± 2) — (8.1 ± 0.4) ■ x. 

5 Conclusions 

We have investigated the hadronic shower longitudinal and lateral profiles 
on the basis of 100 GeV pion beam data at incidence angle = 10° at 
impact points Z in the range from —36 to 20 cm. 

• Some useful formulae for investigating of lateral profiles have been 
derived: 

— the integral expression (^) for radial density $(r) as a function 
of marginal density f{z), 

— the formula (|18|) for radial density $(r) and the formula (16) and 
(p!7|) for cumulative function F{z) in case of three exponential 
form oi f{z) (P). 

• We have obtained for four depths and for overall calorimeter: 

— the energy depositions in towers, E{z); 

— the cumulative functions, F{z)] 

— the marginal transversal densities, f{z); 

— the underling radial energy densities, $(r); 

— the containment of a shower as a function of radius, /(r); 

— the radii of cylinders for given shower containment; 

— the fractions of "electromagnetic" part of a shower; 

— the differential longitudinal energy deposition AE/Ax; 

— the three-dimensional hadronic shower parametrisation. 

We have compared our data with relevant data for conventional iron- 
scintillator calorimeters, SPACAL lead-scintillating fiber calorimeter and 
Monte Carlo calculations. Our longitudinal profile agree with the ones for 
a conventional iron-scintillator calorimeters and Monte Carlo prediction. 
Our lateral profile is not agree with the Monte Carlo prediction. 
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The three-dimensional hadronic shower parametrization for iron-scin- 
tillator calorimeter have been obtained. This parametrisation is important 
in fast Monte-Carlo simulation for ATLAS calorimetry. 

6 Appendix 1. 

Solution of Abelian equation. 

It is stated above that the marginal density distribution f(z) is con- 
nected with the radial energy density $(r) by relation (p!o|) . This integral 
equation can be reduced to the Abelian equation |T^. Here we show how 
to solve the equation (|10|) and to obtain the expression ([Tl|) . Let ^ = 
and 7] = so that the equation ( [T0| ) becomes 



If we multiply ([25| ) on y/rj — a and obtained product integrate over drj 
in limits [a, oo] then 




(25) 




(26) 



Using the following relation 




(27) 



we get 



drj 



(28) 



We should be able to use the fact that 




(29) 



and write 




(30) 
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Let us return to former variables z, r and a = r^. Then and 
modified into 



r \f . 



dz = 71 ^{r)rdr (31) 



Differentiating equation (0) over r we get 



$f ) - ^ ^ f{z)dz'^ 

IT dr'^ . - r^) ' 



7 Appendix 2. 

The effective nuclear interaction length, the effective 
radiation length and the effective Moliere radius. 

In reality our calorimeter represents itself the complex structure of var- 
ious materials and it is necessary to know the effective nuclear interaction 
length the effective radiation length [X^^^) and the effective Moliere 

radius {R% ). 

We calculated these quantities for our calorimeter. For calculating X^ff 
and Xl^^ we used the algorithm suggested in [|I^ so, 

yeff _ ^Fe + Xsc + Xw + Xa ^^^^^ 



Xpel ^Fe + Xsc/ ^Sc + Xw / + Xa/ ^/ 



where xpe, xsc, xw, xa are the volume fractions of the respective mate- 
rials (Fe, scintillator, wrapping, air) in a period of 18 mm thick of the 
calorimeter, Aj are the corresponding interaction length. 

The effective radiation length, X'^J^ ^ was also calculated by the formula 
(|33|) replacing the corresponding Aj values by Xo^i values. The Moliere radii 
is equal [0 

Rm = XqEs/Ec (34) 



and for a mixture of compound 



^ (35) 



Rm Eg ■ Xi 



where the fraction by weight, E^ = 21.2 MeV is the scale energy, 

Ecj are the critical energy. Critical energy for the chemical elements with 
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Table 6: 



Material 


Xj 


Xj, mm 


X^, mm 


Fe 


14/18 


168 


17.6 


Sc 


3/18 


795 


424 


Wrapping 


0.2/18 


Xsc 




Air 


0.8/18 


747000 


304200 



the atomic number of Z is equal Ec = 610 MeV/ [Z + 1.24). The values of 
the Xj, Xj and are given in Table ^ 

In the Table ^ the results of the our calculations are given. In the 
fourth column the corresponding values for SPACAL are also shown 
for comparison. The corresponding values for basic materials of these 
calorimeters (Fe, Pb) from 141 are also given. 



Table 7: 





THECAL 


SPACAL 


eff, mm 


Fe, mm 


eff, mm 


Pb, mm 




22.4 


17.6 


7.2 


5.6 


Rm 


20.5 


16.6 


20. 


16.3 


Xp 


206. 


168. 


210. 


171. 


Xn 


251. 


207. 


244. 


198. 



It turns out that the values A*^-^-^ and R^f/ for these two calorimeters 
are approximately equal. The effective nuclear interaction length for pion 
was also calculated by using relation X-^/Xp (Appendix 3). 



8 Appendix 3. 

The nuclear interaction length for pion. 

Usually in hadronic calorimetry the nuclear interaction length (A/) 
given in Review of Particle Physics [|1^ is used. It is the mean free path 
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for protons between inelastic interactions, calculated using the expression 



A, = A/{NAaip), (36) 

where A is atomic weight, A^^ is Avogadro number, o"/ is the nuclear in- 
elastic cross section, p is a density. For Fe (cr/ = 703 mb) it amount to 168 
mm. But sometimes the pion interaction length is needed as in our case. 
We calculated the pion and proton interaction lengths and compared with 
11^ and with the some experimental data p9|, pSf. For this we used the 



absorption cross section of pion and proton on nuclei in the range 60 280 
GeV/c measured by ||30|. Unfortunately in this work absorption cross sec- 
tion of TT and p on Fe nuclei are not measured. Therefore we used the 
nearest cross section for Cu (see Table Transformation from Cu to Fe 
was performed using the expression 

(TPe = aCuiApe/AcuT- (37) 



Table 8: Measured absorption cross sections and parameter a for Cu 
target. 



E 


60 GeV 


200 GeV 




cr, mb 


a 


cr, mb 


a 


7r+ 


627 ± 19 


0.764 ± 0.01 


629 ± 19 


0.762 ± 0.01 


p+ 


764 ± 23 


0.719 ± 0.01 


774 ± 23 


0.719 ± 0.01 


Avr/ Ap 


1.22 ±0.05 


1.23 ±0.05 




-0.045 ±0.014 


-0.043 ±0.014 



In Table || the measured cross sections for tt+Cm interactions at 60 and 
200 GeV together with parameter a are given. Using expressions (|36| ) and 
(|37|) cross sections and nuclei interaction lengths are presented in Table ^ 

As can be seen from Table ^ the Fe cross sections and nuclear inter- 
action lengths do not depend from energy within errors at 60 -v- 200. The 
mean over energy range values are A,r = 207 ± 7 mm, Ap = 168 ± 5 mm, 
Ajr/Ap = 1.22 ±0.05. The value of Ap coincides with that given in |14 
value of Att may be compared with the measured in p5[ value A^r 
mm and in |^ value A^ = 190 mm. 



The 
200±10 
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Table 9: Calculated cross sections and nuclei interaction lengths for Fe 
target. 



E 


60 GeV 


200 GeV 


a-jr, mb 


568 ± 18 


570 ± 18 


(Tp, mb 


696 ± 21 


705 ± 23 


A^, mm 


207 ± 7 


207 ± 7 


Ap, mm 


169 ± 5 


167 ± 5 


Att/ Ap 


1.22 ± 0.05 


1.24 ± 0.05 



The value measured in presumably have be corrected on iron equiv- 
alent length (including scintillators) and then amounts to 210 ±10 mm and 
better agrees with ours. 

Using the expressions ( p6D and ( p?!) we obtained the ^-dependence of 
the ratio 

AJAp(A) = A^/Ap(Cn)(AMc„)"''-"^ (38) 

where Op — = 0.045 ± 0.014 (Table P). For Pb nucleus it corresponds to 
\/\p{Ph) = 1.16 ± 0.05. Therefore, it is mistaken to use the value of the 
ratio A^/Ap(P6) = 1.5 as have been made in when analysing of data 
from lead scintillating fiber spaghetti calorimeter. 
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